Dissipative solitons which cannot be trapped 
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In this paper we study the behavior of dissipative sohtons in systems with high order nonlinear 
dissipation and show how they cannot survive under the effect of trapping potentials both of rigid 
wall type or asymptotically increasing ones. This provides an striking example of a soliton which 
cannot be trapped and only survives to the action of a weak potential. 
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Introduction.- Solitons are self- localized nonlinear 
waves which are sustained by an equilibrium between 
dispersion and nonlinearity and appear in a great va- 
riety of physical contexts [ij. Solitonlike localized states 
also arise in dissipative systems driven far from thermal 
equilibrium such as Jiydrodynamics 0], granular media 
nonlinear optics and other fields 
These structures are re- 
ferred to as "dissipative solitons" and are sustained be- 
cause of an interplay between dispersion, nonlinearity, 
gain and losses. 

A new family of multidimensional dissipative (bright) 
solitons has been recently observed and studied in the 
context of the propagation of optical beams in media with 
multiphoton absorption [la, llll llM ll^ ■ These solitons 
have slowly decaying tails and have a nontrivial flux of 
energy from their tails to the soliton center where energy 
is absorbed leading to a stable nonlinear object. In some 
sense they represent a nontrivial extension of the well- 
known non-diffracting optical beams j2ni |. 

These dissipative soliton-like structures have been also 
studied in different contexts such as three-dimensional 
scenarios j2ll and should also exist in other physical 
scenarios ruled by analogous physical phenomena (and 
model equations) |2^ . 

In this letter we will describe a highly counterintu- 
itive property of these solitons which is that they can- 
not be trapped by the action of an attractive potential, 
an affirmation which we will make more precise later. 
This fact is very striking since attractive potentials use 
to enhance the stability of solitons. In fact, it is well 
known from Quantum Mechanics that attractive poten- 
tials alone without any help of nonlinear interactions are 
able to support stationary solutions. Even unstable non- 
linear structures such as the Townes soliton 
bilized by an external confining potential [2. 

Mathematical model.- As a mathematical model of soli- 




tons under high order dissipation we will use the nonlin- 
ear Schrodinger equation (NLS) 

dib 1 

t-^^--A^ + V^P + g\,p\'^ + tj\^\^i:, (1) 

on R^, where g and 7 > are real parameters (non- 
linearity and dissipation coefficients) and is a real 
function describing an external potential acting on the 
system. Eq. ^ models the propagation of optical 
beams in media with multiphoton absorption [THIT?, 21j . 
the mean field dynamics of Bose-Einstein condensates 
with losses associated to inelastic three-body collisions 
Q2i,Mi,M,28^2^-,M,ILM and other physical phenom- 
ena [23 . The mathematical studies on this equation have 
focused on the study of the concentration mechanisms 
leading to the so called "super-strong" collapse phe- 
nomenon (which however, does not involve the formation 
of a singularity in the mathematical sense) I2H (sS, ^S^l • 

Stationary solutions.- Stationary solutions of Eq. 
with 7 = 0, and = have been discussed extensively 
since they correspond to the simplest model of nonlinear 
optical beams and classical superfluids. When 7^0 they 
have been studied in Refs. llTl I21I, l2a| • Let us write 



ip{r, t) = A(r, t) exp[i0(r, i)], 
with A{r, t) > 0. Then Eqs. Q become 



dt 



(2) 



(3a) 



dt 



l^A-{\/<t>f 



gA^ - V. (3b) 



Stationary solutions of Eq. satisfy dtA = 0, and 
(f>{r,t) = ipir)-Xt. 

It is obvious that N = ||-0||2 = Jjg.2 IV'P = /112 A'^, can- 
not change over stationary solutions. However, let us as- 
sume that N is finite for stationary solutions, then using 
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Eq. (P) we get 



dN 
~dt 



= -27 / < 0, 



(4) 



which is a contradiction. This impUes that N cannot 
be finite for nontrivial stationary solutions of Eq. 1^. 
This fact is also known to be characteristic of the linear 
non-diffracting optical beams [20|. 

The existence of a nontrivial phase structure ^p{r) is 
a feature of the case 7 7^ which is not present in the 
conservative NLS. From Eqs. I|3a() and (|3b|l 
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'2.1 A^, 
\A-gA^. 



(5a) 
(5b) 



In this paper, in order to deal with physically meaning- 
ful solutions, we will study solutions of Eqs. ^ which 
decay asymptotically for large values of r = ||r||, i.e. 
A(r), ^p{r) when r — > 00. 

When the amplitude and phase are radially symmetric 
functions A{r) — R{r), ip{r) = <I>(r), Eq. can be 

integrated by using the divergence theorem and we get 



Introducing JB)) into Ij5b|l we can write 



(6) 



R" - -R' + 2VR 



2XR - 2gR^ - 



47^ 



pRipfdp 



(7) 



complemented with the boundary conditions R{r) — + as 
r — > 00 and i?'(0) — 0. Eq. Q is a nonlinear eigenvalue 
problem with a nonlocal term. 

Let us look for solutions satisfying that 



g = 27/ ri?(p)«dr = J||i?||ie(g.), (8) 



is finite (mathematically, this means that R S L^{M.'^)). 
In this case Eq. ® allows us to obtain $'(r) « 
Q/ [ri?2(r)] , for r > 1, and 



-R" 



-R' + 2VR = 2\R - 2gR^ - 



R^r^' 



(9) 



System without external potentials (V=0).- When A > 
and V — Q it was shown in Ref. 0| how for large 
r values the leading order approximation to the solu- 
tion is R{r) ~ When A < there are no sta- 
tionary solutions (this was conjectured in Ref. ^^'^ 
will be proven here for later convenience. First we write 
A{r)e'^f^^'> = u + iv, then Eqs. 1^ become 



-Av = Xv - g{u'^ + v^)v - j{u'^ + v^)'^u. (lOh) 



Let us now consider Eqs. H10|l on a domain given by a 
generic ball Bp with boundary conditions u = v = on 
dfl and u, v to be radially symmetric solutions. Multi- 
plying the first of these equations by u, the second by v, 
integrating in Q and applying Green formulae we get 



np (uUr 



X 



- V 



Bo 



9 1 (w' 



,,2\2 



(11) 



Since R > and i? ^ at infinity, we can always take 
a sequence pk 00 such that Rr{pk) < 0. However 
{R'^)r — 2{uur + vvr), therefore the l.h.s. of Eq. l(TT|) 
is positive for any Bp^ . As a conclusion, if g > there 
are no radially symmetric solutions. 

When A < 0, (? < and V = Q there are not radially 
symmetric solutions with amplitude in L^(S?). By com- 
parison results on the Eq. ® , it is not difficult to prove 
that, when r ^ 1, R{r) > C exp(\/— A r) for some C € M, 
which is unbounded when r 00 provided A < 0. 

In summary, when A < there are neither radially 
symmetric solutions at all for g > nor solutions sat- 
isfying the finiteness requirement for Q (mathematically 
A G L^(R2) norm) for g < 0. Thus, we can ensure that 
there are no physically relevant stationary solutions when 
A < 0. 

Nonexistence of dissipative solitons in systems with 
rigid walls.- Now we move to the main point of this pa- 
per, which is showing that dissipative solitons of Eq. 
cannot be trapped by strongly confining potentials, or 
equivalently proving that stationary solutions of Eq. 
do not exist for those potentials. 

To do so we first consider the case of an infinite po- 
tential barrier placed on the boundaries of a finite region 
n. Energy conservation then imposes the condition that 
u{r) — for r e dil. 

We start from Eqs. ifTUIl . multiply Eq. HlUa|l by v, Eq. 
(jlObp by u, and integrate over to get 

- / Vu • = A / uv — g f [u^ + v'^)uv 
2 Jsi Jn Jn 

+7 / {u^ + v'^fv'^, (12a) 

— ( Vm • Vv — X f uv — g f {u^ + v'^)uv 
2 Jn Jn Jn 

,2 , „,2^2„,2 (^2b) 

The difference of Eqs. (|12a|) and (jl2b|l leads to 



-7 / {u' + v^u' 
Jn 



7 / [u'^ + v^f^O, 
In 



(13) 



from where we conclude that either 7 = or the only 
solution is the trivial one u = v = 0. Thus dissipative 
solitons cannot exist when confined to a finite region. 
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We can provide some extra results valid for more gen- 
eral geometries of the domain fi. Let us consider again 
Eqs. 10 . Multiplying Eq. ifTUajl by u, Eq. l(ToH)l by v, 
integrating in Q and applying Green formulae we get 

2 Jo. Jn Jq Jn 



In Jn Jn Jn 

Adding both equations we obtain 

1 / \Vu\^ + \\7v\^ = X [ u'+v^^g [ {u^ + vY- (14) 

2 Jn Jn Jn 

Poincare's inequality on bounded domains Q says that if 
di is the first eigenvalue of the Laplacian operator on 
with Dirichlct homogeneous boundary conditions, then 



Si = inf 



In IV^I^ 



(15) 



This inequality also holds for unbounded domains Q, 
which are bounded along one of the directions, which 
allows us to extend our results to systems with strong 
confinement along a single direction. In any of those sit- 
uations we will have 

^ (u^ + v^) <l [ |VuP + |V«P 



X {u^ + v')^g {u^+v^r, (16) 
Jn Jn 



from where we find that 



<-9liu'+v'r, (17) 
n 



so that we get that when g > the negativity of the 
second term implies that A > 5i/2, thus if A < Si/2 
the only solution is u = w = 0. Moreover Si{fl) 
when \ — !■ 00, which matches our previous results on 
radially symmetric stationary solutions. However, here 
the result is proven rigorously for any type of solution 
(not only symmetric ones) in hard-wall type potentials 
corresponding to semi- infinite regions. 

Systems with "strong" confinement .- We have already 
shown that dissipative solitons of Eq. ^ cannot be con- 
fined to finite spatial regions by hard-wall type poten- 
tials. In what follows we will study the effect of nonsin- 
gular potentials. To do so we will first start by considering 
unbounded asymptotically increasing potentials i.e. po- 
tentials with y ^ 00 as r — !■ 00, such as those of the 
form < V{r) ^ with p > 0. This is the case of appli- 
cations to Bose-Einstein condensates in magnetic traps, 



strong confinement and in the linear (5 = 0, 7 = 0) case 
have an infinite number of bound states. 

Let us then consider the version of Eq. (|11|) including 
the potential term. The same type of arguments used for 
the case y = 0, A < can be used to prove that if 5 > 



X > 



IbR' 



(18) 



As stated before, stationary solutions must have neces- 
sarily Jg i?^ ^ 00 when p ^ 00. For asymptotically 
increasing potentials there is some r = such that for 
all r > r*, y(r*) > A and then both the numerator and 
denominator of Eq. (|18|l diverge and using L'Hopital rule 



lim 



IbZ^ 

IbR' 



lim 



i?2 



(19) 



therefore A > limp^oo Vip) = 00 ^^nd as a conclusion, 
there are no radially symmetric solutions. 

When (? < 0, there are no radially symmetric solutions 
with amplitude in L^{R'^). The reason is again that, when 
r ^ 1, R{r) > C exp(V— A r) for some constant C, which 
is unbounded when r — s- cx). This excludes the possibility 
of existence of physically meaningful stationary solutions 
for both signs of the nonlinear coefficient g and implies 
that dissipative solitons of Eq. cannot exist when an 
asymptotically increasing potential is present. 

Systems with weak confinement.- We will now see how 
dissipative solitons survive when a bounded potential 
(weak confinement) is present, i.e. V{r) —^ Voo when 
r ^ 00. In this situation, arguing as in the previous 
case. 



A > inf 



IbZI 

IbJ^ 



> lim inf V 



(20) 



where usually V(r) 



These potentials provide a 



thus implying the possibility of existence of solitary waves 
only for values of the eigenvalue A satisfying A > 14o • In 
fact, the asymptotic analysis of the problem for r 3> 1 in 
this case is the same presented in Refs. 0, and we 
can expect the existence of solutions which will be close 
to those of the case V = but with local deformations 
in the region where V has appreciable values. 

Physical implications and soliton mobility.- These re- 
sults have practical implications for the experimental 
generation of dissipative solitons in Bose-Einstein con- 
densates with appreciable losses due to three-body col- 
lisions. In realistic condensates, due to the condition 
< cx) wc can only construct dissipative structures 
close to the stationary solutions but with a finite lifetime 
|22l |. In fact, numerical simulations of blow up phenom- 
ena with high order losses and finite number of particles 
show the formation of multi-ring structures whose shape 
resemble that of dissipative solitons [2^ E^. However, 
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our results imply that the magnetic confining potentials 
of the form V{r) ^ are not appropriate to confine the 
atomic cloud if these dissipative solitons are to be ob- 
tained. In this sense it would be more convenient to use 
a laser beam to generate an optical dipole type trap, of 
the form V{r) ^ e'~^ thus being compatible with the 
existence of stationary structures. 

Could a localized potential be used to control these 
solitons? The problem of soliton mobility is an inter- 
esting one in many different contexts |3{tI | . In our case 
an optical dipole trap (or any other "localized poten- 
tial") cannot probably be used to move the soliton. The 
reason is that these dissipative solitons arise as a global 
equilibrium between an inward fiux of particles due to 
the nontrivial phase structure and the central peak at 
which particles are dissipated So, even if a local 

potentials could be used to move the central peak, the 
refillin g rn echanism would lead to a new soliton at the 
center Il8| while the moving part would experience a fast 
decay due to losses. The only possibility for moving this 
soliton would be to use a bounded potential (to ensure 
existence) having a very slow decay ensuring a global, 
although small, action on the dissipative soliton (i.e. on 
the same scale of the soliton). An example could be a 
potential of Yukawa type V{r) ~ [1 — exp(— r)] /r. 

Conclusions.- In summary, we have studied a class dis- 
sipative solitons in models which have received a lot of 
attention in the last years in differe nt phys ical contexts 

[IlIilillllillSllllMEilMEllIi 

Our main result is that contrary to the intuition and 
the observed typical behavior of solitons these dissipative 
solitons cannot be trapped by using strongly confining 
potentials. This fact represents a notable exception to 
the typical behavior of solitons in the context of systems 
ruled by nonlinear Schrodinger equations. We have also 
discussed some physical implications of our findings and 
made some conjectures on the mobility of these solitons 
on the basis of our findings. 
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